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We present an operator-algebraic approach to deriving the low-lying quasi-degenerate spectrum 
of weakly interacting trapped N bosons with total angular momentum hL for the case of L/N <C 1, 
demonstrating that the lowest-lying excitation spectrum is given by 27 gn?,(nz — l)/34, where g is 
the strength of the repulsive contact interaction and 713 the number of excited octupole quanta. Our 
method provides constraints for these quasi-degenerate many-body states and gives higher excitation 
energies that depend linearly on N. 
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I. INTRODUCTION 



Realization of Bose-Einstcin condensates (BECs) in 
atomic gases has stimulated explosively growing in- 
terest in several subfields of physics. With unprecedented 
controllability in creating, manipulating, and probing 
systems, many previously "gedanken" experiments have 
become or are expected to become a reality [|J . Another 
characteristic feature of gaseous BEC systems is their 
good isolation from the surrounding environment, as the 
system is suspended in a vacuum chamber by a magnetic 
or an optical trapping potential. When the potential is 
axisymmetric, the projected angular momentum (AM) 
on the symmetry axis is conserved and can take on any 
integral multiple of fi. In contrast, in bulk superfluid sys- 
tems, the container plays the role of a highly dissipative 
environment so that the system can sustain AM only at 
thcrmodynamically stable values, resulting in the quan- 
tization of circulation. In this respect, gaseous BEC sys- 
tems are more like nuclear systems than bulk superfluid 
systems, and it is of interest to study the energy spec- 
trum of a weakly interacting BEC system as a function 
of AM, namely the yrast spectrum 

The lowest-energy state of an isolated system subject 
to a given AM, fiL, is called the yrast state 0. When 
the system is in the yrast state, all kinetic energy is used 
for rigid-body rotation, so the system is at zero tem- 
perature. This allows one to measure fine energy levels 
— the yrast spectrum — close to the ground state. Re- 
cently, the yrast state has been extensively discussed in 
the context of gaseous BEC systems (5]j|,|l^ |l6| . Mottel- 
son has pointed out that for repulsive contact interactions 
the yrast spectrum is dominated by quasi-degenerate 
quadrupole and octupole excitations ||], and Kavoulakis 
et al have shown that the interaction between octupole 
excitations makes the excitation energy of this mode 
slightly higher than that of the quadrupole mode JTo[ . 
Bertsch and Papenbrock jl2j have found numerically that 
for < L < N, the energy of the yrast state is given by 



E = LTico + gN(N — 1 — L/2), where oj is the frequency 
of the confining potential, g the strength of the interac- 
tion, and N the number of bosons, which we shall assume 
to be constant. The corresponding eigenstate has been 
shown to exist in Refs. p3|-p^|. Recently, we have found 
numerically that the lowest-lying yrast spectrum arises 
from the pairwise repulsive interaction between octupole 
excitations and that the corresponding interaction energy 
is given by 1.59 17713(713 — l)/2, where 773 is the number 
of octupole excitations |lq| . In this paper, we present an 
operator-algebraic approach ]l7| to analytically deriving 
these results in a systematic way, demonstrating that the 
interaction energy between octupole excitations is given 
by 273773(773 — l)/34, where the coefficient 27/34 agrees 
precisely with the numerically obtained value [fL6| . Our 
method also provides constraints for the quasi-degenerate 
many-body states and gives higher excitation energies 
that depend linearly on N. 

This paper is organized as follows. Section || describes 
the model u nder consideration and formulates the prob- 
lem. Section III describes collective modes and discusses 



why they behave as well-defined quasiparticles when the 
system undergoes the Bose-Einstein condensation. Sec- 
tion ^ derives the yrast spectrum based on our operator- 
algebraic approach, and Sec. ^ concludes this paper. 



II. FORMULATION OF THE PROBLEM 



We consider a two-dimensional system of harmonically 
trapped N bosons interacting via a contact <5-function in- 
teraction. The Hamiltonian of the system is given by 



H = 
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where M is the mass of the boson, lu the frequency of the 
confining potential, and g the strength of the interaction. 
In the second-quantized form, Eq. ([I]) is written as 



H 



d l z 



¥{z)h(z)i>{z) 



2nhg 
Muj 



where the prefactor is determined so as to satisfy the nor- 
malization condition J d 2 z \(f) (z)\ 2 = 1. The wavefunc- 
tion 4>m(z) for the state with magnetic quantum number 
m is given by 



&(z)&(z)&(z)${z 



where z = 



iy and 
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denotes the single-particle Hamiltonian. 

Throughout this paper we consider the case of very 
weak interactions, i.e., \g\N <C Tito. It is then convenient 
to take the eigenfunctions of h(z) as a basis set for ex- 
pansion of the field operator. We introduce two sets of 
creation and annihilation operators, 



z 

2d 

z 
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dz' 
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dz* 
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(4) 
(5) 



where d = y/Ti/Muo is a characteristic length scale of the 
ground-state wave function of h(z), and a + (or d_) anni- 
hilates a quantum with magnetic quantum number 1 (or 
-1) and satisfies the canonical commutation relations 

[«+)«+] = [a-,a f _] = 1, [a+,a_] = [d+,d t _] = 0. (6) 



In terms of these operators, h is diagonalized as 



a d_ 



!)• 



(7) 



The vacuum state |0) is defined by a+|0) = d_|0) = 0. 

The operator corresponding to the magnetic quantum 
number m is given by d^a + — a_a,-. Because of the 
symmetry between d + and d_ , we may consider only the 
case of m > with no loss of generality. Moreover, when 
Two ^> \g\N, we may consider only those states {IvP)} that 
satisfy 



•l*> = 



(8) 



because for a given total AM, states that do not satisfy 
Eq. (^J) are at least 2Tno higher in energy than those sat- 
isfying the equation. Physically, Eq. (||) suggests that 
no particle rotates in a clockwise direction. In what fol- 
lows, we shall restrict ourselves to the Hilbert subspace 
constrained by Eq. (||), which corresponds to the lowest- 
Landau-levcl approximation in the theory of the frac- 
tional quantum Hall effect. 

Within this approximation, a complete set of basis 
functions can be constructed as follows. Using (z|a+|0) = 
and Eq. (@), we obtain 



>(*) 



-1°) 



1 
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/7rm! d m+1 
and satisfies the orthonormal conditions 



(z\0) 
(10) 



d 2 z (j)* m {z)(f> n (z) = 8„ 



(11) 



It follows from Eq. (jlTj) that \z) is a coherent state de- 
fined by 

\z) = (V^d)- 1 exp[{z*al - za+)/d] |0). (12) 
Expanding the field operator in terms of 4> m {z) as 

*(z)=5^6 m ^ m («), (13) 

m 

and substituting this into Eq. (|2|), we obtain 
H = TwoL + g ]T V m (14) 



mi,- ■ ,7714 



where 



(15) 



gives the total magnetic quantum number and 
is the matrix element of the interaction given 



as |l 



V, 



J 7ii i +m2 ,m3+m4 



(mi + m 2 )\ 



2 m i+™2 y/mi\m 2 \ m 3 !TO 4 ! 



(16) 



In Eq. ( |lj| ) the term TitoN is dropped because the total 
number of bosons 



(17) 



is conserved in our system. 

Because L is also conserved in our system, we shall 
henceforth focus on the second term on the right-hand 
side (rhs) of Eq. (Q), which will be denoted as V = gU. 
Our primary goal is to diagonalize U in operator form 
within the Hilbert subspace of non-negative m and un- 
der the two constraints 



L, N = N. 



(18) 
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III. COLLECTIVE MODES 



Mottelson has pointed out that low-lying excitations 
from the yrast line can be described in terms of collec- 
tive modes excited by (d~ x /y/N\\) Xw=i z j' where A is 
a positive integer. In the second-quantized form, the cre- 
ation operators for the collective modes are written as 



Q{ = / d 2 z&(z) 



\J \+ m C\ b 



z x ^(z) 



£ V x+mCx b m+x b m , 

m—0 



(19) 



where m C n = m\/n\(m — n)\. The commutation rela- 
tion between Q\ and the interaction Hamiltonian V is 
calculated to be 



[V,Q> 



2.9 



(mi + m 2 )! 

'ATA! ^ 2 m i+ m 2+ A Jm\\rn-x\mz\m±\ 

TTh\ , • ■ ■ . TYI4. 



E 



1114! 



(mi + m 2 + A)! 
(m 4 — A)! (mi + m 2 )! 



(20) 



It follows from this relation that except for A = 1, Q\ 
does not commute with V. The collective mode excited 
by Q\ (A > 1) is therefore not, in general, a well-defined 
quasiparticle. The only exception is the dipole operator 
Qi, for which the rhs of Eq. ( po| ) vanishes. The dipole 
mode is therefore a well-defined collective mode, and for 
any L, an eigenstate of the Hamiltonian (|l4| ) may be 
constructed as j|D 

\^) = {Q\f\L = 0), (21) 

where \L = 0) denotes the exact ground state of the 
Hamiltonian ( fl4| ) with zero total AM, and hence V\L = 
0) = gN(N — 1)| L = 0). The corresponding eigenenergy 
is given by 

^dipole 



TiuL + gN{N - 1) 



(22) 



because 



H \^ L ) = (hcuL + V) (Q\) \ L = 0) 

= huoL |* L ) + (Q\) L V\L = 0) = E dipolc 

The state |\£ ) is the lowest-energy state of the Hamilto- 
nian with AM L for the case of an attractive interaction 
(g < 0) and the highest-energy state for the case of a 
repulsive interaction (<? > 0). 

It should be noted that while Q\ and commute for 
any A and fi, Q\ and CK do not. In fact, 



m+A m+A 



m=0 

— \J m+\C\ 



(23) 



where 9(x) is the Heaviside step function, which takes on 
0, 1/2 and 1 for x < 0, x — and x > 0, respectively. 
The commutation relation (|2^) implies that the collec- 
tive modes carrying different angular momenta do not, 
in general, behave independently. For pb = A, Eq. fl23| ) 
reduces to 



h,Q{] 




for A = 1; 
for A = 2; 

for A = 3, 



(24) 



where constraints ( |l8[ ) arc used in deriving the second 
equality. The commutator [Q\, Q\] in Eq. ( p4| ) gives con- 
stants for A = 1,2, indicating that the dipole (A = 1) 
and quadrupole (A = 2) modes behave like bosons. In 
contrast, the commutator for the octupole (A = 3) mode 
does not give a constant, indicating the anharmonicity of 
this mode. This anharmonicity leads to the interaction 
between octupole excitations, as shown in Sec. 



IV 



The energy minimization of V must be performed sub- 
ject to constraints (|l8|). In the case of L/N <C 1, which 
we shall consider below, these constraints require that 



n 



0{N), n m =O(N ) (m>l), 



(25) 



where n m = b\j3 m and 0(N k ) denotes terms of the order 
of or less than N k . In this case, Q\ may be simplified as 



N 



blbx, 



(26) 



and the commutation relations in Eq. (|23]) reduce to 



[Q\,Qi]=5x» + G(-). 



(27) 



The collective modes, therefore, behave as almost inde- 
pendent bosons for L/N <C 1, with small inter- mode in- 
teractions on the order of L/N. These features are re- 
flected explicitly in the yrast spectrum, as shown in the 
next section. 



IV. YRAST SPECTRUM 



Our basic idea for the derivation of the yrast spectrum 
is to rewrite the interaction Hamiltonian V = gU under 
constraints ( |l8[ ) so that the excitation spectrum mani- 
fests itself. We first note that U can be rewritten as 



U — £ Mn-h> 



(28) 



m=0 
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where 



Substituting these expressions into Eq. (E8h yields 



A m — 2 2 2_j WmPk bkb m ~k 



(29) 



k=0 



may be regarded as the annihilation operator of a boson 
pair with the total magnetic quantum number m. As we 
shall consider only the case of g > 0, the energy mini- 
mization for U may be exploited to estimate the order of 
magnitude for each term. 

The conditions (Eq) allow us to expand A m as 



H b m + O(N°) (m> 1). (30) 



Substituting these expansions into Eq. 



we obtain 



U = N(N - 1) - 27V V (1 



m— 1 



l )n m + 0{N 1 ' 2 ) (31) 



where the constraint N — N is used to eliminate no- Be- 
cause h 2 and h 3 appear on the rhs of Eq. ( pl| ) in a form 
proportional to 2fi 2 + 3n-3, the constraint L = L may be 
used to eliminate both of these terms simultaneously. We 
thus obtain 



U = N IN- 1 



L 



+N 



- + E(y- 2 + 2 



m=4 



©(TV 1 / 2 ). 
(32) 



The minimum of U in Eq. (|32j) is attained when 

hi =(D(N- 1/2 ), n m = (D(N- 1/2 ) (m>4). (33) 

These conditions are compatible with the constraint N — 
N because no = 0(N°), and with the other constraint 
L = L because it can be satisfied by an appropriate choice 
of n 2 and n 3 . We note that the first term in Eq. ( |32"|) gives 
the yrast line Jl2] |. 

Because the yrast line is dominated by the m = 2 and 
m = 3 modes, the low-lying excitations from the line due 
to binary interactions ( |l4| ) should involve modes with at 
least up to m = 6. We therefore use A m up to m = 6 
without approximation: 

p 

A =bl, Ax=V2b b u A 2 = b b 2 + -=, 

v2 



^3 = 

A 5 = 
A 6 = 



bpb 3 
V2 
b h 
V8 
bpb 6 
4 



A 4 = 



b b 4 



bib 3 



la. 



5 £ t ? 
-01O4 + -^-o 2 b 3 , 



/ 3 £ i 

■xhb 5 + 



b 2 b 4 + ^bj. 



(34) 



L \ 27 

U = N [ N - 1 - -\ + —n 3 (h 3 - 1) + R + S 



N . 



TV 



5iV 



._B[B l + -B[B i + —Bl 3 B b 



17 N 



B|5 6 . (35) 



Here R is given by 

E (~ 7^ ' ' + 

OO OO 



E 

fc=6 



fc - 3 



-nin k 



E 

fc=5 



fc-2 



-n 2 n fe 



A: + I — 4 „ „ 5.„ 
2 , ^ n fe n; - 2(n 3 + l)n 4 - -(n 4 + l)n 5 

A:=3i=fc+1 

-i(^Sp|6 2 S 5 + 2V5bfb 3 b 5 + h.c.) 
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[23riin5 + 32h 2 n 4 



~(uVlOb{blb 2 b 4 - VWblblbj + l2V3blblb 2 3 + h.c.)], 

(36) 

where h.c. denotes Hermitian conjugates of the preceding 
terms, S is given by 



S = y~] ( A\ n A m + m 4 n n m , 



m— 7 



(37) 



and B m (m = 1, 4, 5, 6) are defined by 



B x = 



1 

T^E^I-i^' 
fc=l 



6 6 4 - 26i6 3 



N 



b\ , 



Br, 



Ba 



N 



6 6 5 ^6164 + \ -b 2 b. 



N 



bobe 



^/5 

I x ^ A 



17 - 



(38) 



Note that no approximation is made in obtaining 
Eq. (H). It follows from (H) that R = 0(N-^ 2 ). The 
operator 5* is positive semidefinite. To minimize it, we 
consider the states that satisfy 



ft fc |*) =o(N-V 2 ) (fc>7), 



(39) 



so as to give = o(iV ), where o(iV Q ) denotes 

terms whose order of magnitude is less than N a . Under 
the same condition we obtain 

U = N(N-l--\+-n 3 (n 3 -l) 



N . 



5iV . 



-jB{B 1 + T B^ + —B^B s 
17 N 



16 



BlB 6 +o(N°), 



(40) 
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which constitutes the primary result of our paper. 

Let us define the quasi- degenerate states, {I'l')}, as the 
ones satisfying both condition (B9J) and 



B m \9) =o(AT3) ( m = 1,4,5,6). 



(41) 



The following constraints on these states are then ob- 
tained: 



(B} n B m ) = o(N- r ) (to = 1,4,5,6), 
{n m )=o(N- 1 ) (to>7). 



(42) 
(43) 



The expectation values of U for these states are therefore 
given by 



(U) = N (N - 1 - 



L 



)+§ n 3 (n 3 -l) + o(N°), (44) 



where n 3 is the number of octupole excitations, as 
Q[Qx = n x + OiN- 1 / 2 ) for A > 1. Since the popu- 
lation of the octupole mode with n 3 > 2 raises the en- 
ergy, while that of the quadrupole mode does not Jl0|,|l6| , 
the ground-state configuration is given by ri2 — L/2 and 
n 3 = when L is even and by n 2 = (L — 3) /2 and 713 = 1 
when L is odd. These results are valid up to the order of 
TV , as Eq. (^J) is valid up to the same order. 

As seen in Eq. (^), the quasi-degenerate octupole 
spectrum lies above the yrast states. Although the oc- 
tupole mode couples with other modes through 



07 

[b 3 ,U}--btbt 



-b 2 B 1 



A7J|+ - ^-fb\B 4 



5 b 2 B 5 + 



-btB e , 



(45) 



the expectation value of the rhs over the quasi-degenerate 
state satisfying Eq. (|4j) is on the order of o(N°). The 
quasi-degenerate spectrum 27n3(ri3 — l)/34 is therefore 
valid up to the same order. 

The representation of U in Eq. ( f40| ) not only describes 
the quasi-degenerate yrast spectra, as discussed above, 
but also explains higher excitations whose energies are 
linear in N. Since 

[B m , B n ] = OiN- 1 / 2 ) and [B m ,B& = S mn + OiN^ 2 ), 

(46) 

each B m (to = 1, 4, 5, 6) behaves as a boson-like operator 
and satisfies 



BlB m ^n m + 0{N- 1 / 2 ). 



(47) 



Because of the definitions (|38|) of B m and from (|33J), the 
excitation energies for B±, B5, and Bq are expected to be 



e + e 4 - 2e 2 = gN/A, 

eo + £5 - £2 - £3 = 5gN/8, and 

eo + e 6 - 2e 3 = 17gN/16, 



(48) 



respectively, and the excitation energies for B\ are ex- 
pected to be 

ei + e 2 - e - £3 = gN/2 and 2ei - £ - £2 = gN (49) 

for n\ = 1 and n\ = 2, respectively. Here 

e A = -25^(1 - 2'- x ) (50) 

gives the energy of a single excitation of Q\ (A > 1) j|] 
and eo = 0. In fact, these excitation energies appear as 
the coefficients for B^Bm in Eq. (|40|). 

Here we discuss the physical meaning of Eq. (^Tj) . In 
the BEC the m = state is macroscopically occupied. 
When the BEC has angular momentum L(1<L< N), 
the occupation numbers of m — 2 and to = 3 states 
become of the order of L. Because of the interactions, 
the following scattering processes, whose energy contri- 
butions are of the order of 0(N°), become important: 
+ 2^1 + 1, + 3^1 + 2, + 4^2 + 2, + 5^2 + 3, 
and 0+6 «-> 3+3. The corresponding terms in the interac- 
tion Hamiltonian can be diagonalized by the introduction 
of the quasiparticle operators B m (to = 1, 4, 5, 6) defined 
in Eq. ( pq ) . This procedure is nothing but the renormal- 
ization of the bare operators b m due to the interaction. 
The true vacuum is therefore not the bare one defined by 
b m \^) — (to = 1,4,5,6) but that of the quasiparticles 
defined by B m \^) = 0. This is the physical meaning of 
the constraint (f4l|). 

On the other hand, the contact interaction dictates 
that neither quadrupole (to = 2) nor octupole (to = 3) 
mode can make O(N ) energy contributions together 
with higher angular momentum states (to = 7 or higher) . 
Hence the supplementary constraint ([39|) for to > 7 
should be imposed for the complete definition of the (in- 
teracting) vacuum state. In fact, from Eq. (|3^) b m exci- 
tations are shown to need energy costs which are linear in 
N. We note that because of this renormalization proce- 
dure of the quasiparticle operators, our excitation spec- 
trum achieves the accuracy of o(N°) for quasi-degenerate 
states. Moreover it should be noted that the other low- 
lying excitations can also be explained naturally in terms 
of B m and b m excitations. 

As an independent check of the above analytical re- 
sults, we have numerically performed an exact diago- 
nalization of the Hamiltonian, confirming that the con- 
straints ([l3|) on the quasi-degenerate states are satis- 
fied. The excitations with energies linear in N are also 
found to exist and are shown in Fig. [l] as open circles. 
Here the solid line indicates the yrast line —gNL/2, and 
the dashed, dotted, dash-dotted, and dash-double-dotted 
lines correspond, respectively, to the gN/4, gN/2, 5gN/8, 
and 17giV/16 excitations that are predicted by Eq. (|4C|). 
The broken curve connects the points of collective excita- 
tions Q{\0). We see that the excitations shown in Fig. |l| 
can be well explained in terms of single excitations of 
B m (m = 1,4,5,6). When these higher modes are ex- 
cited, the quasi-degenerate conditions (El]) are no longer 
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met, so that the coupling of these modes to the octupole 
mode through Eq. ( fl5|) becomes significant. This cou- 
pling causes O(A r0 )-order corrections to the excitation 
energies, which are on the order of <D(N). The examples 
of our numerical calculations demonstrating these small 
corrections are shown in Table ||. Since the relative order 
of magnitude of these corrections is on the order of 1/N, 
they cannot be discerned in Fig. [I]. 

Our method provides some constraints on the quasi- 
degenerate states, including the yrast state, through j39| ) 
and (J4l|) . From these constraints, we may gain some in- 
sight into quantum fluctuations occurring in the quasi- 
degenerate states {|^)}. For example, from the condition 
= o(N-i) and @, we obtain 



Repulsive 10000 Bosons 



(blh + V2b\b 2 + v^SJ&a) I*) = o{N°). (51) 

This equation strongly suggests that the interconversions 
between the states with m — 0, 1, 2, 3 (and the intercon- 
versions between quadrupole and octupole modes |16|) 
play an important role in the quantum fluctuations of 
the quasi-degenerate states. We stress here that the rhs 
of Eq. @ is not O{N ) but o(N°). Similarly, from the 
condition B$ | ty) = o(7V - 2), we obtain 



boh 



b 2 b 3 = o(JV°). 
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2 4 6 8 

Total Angular Momentum 

FIG. 1. Energies (in units of gN) of the yrast states 
and those of low-lying excited states with N = 10000 for 
1 < L < 8. Open circles show our numerical results, the 
solid line indicates the yrast line —gNL/2, and the broken 
curve connects the collective modes excited by Q\. The 
dashed, dotted, dash-dotted, and dash-double-dotted lines 
correspond, respectively, to the ffiV/ 4, g N/2, 5gN/8, and 
17gN/16 excitations predicted by Eq. Q. 



It follows from these two constraints that on a mean- 
field level we have (h\) ~ 15(ns)/2. From the con- 
straints B m \^) = o(N~i) (m = 4,5,6), we may also 
have (hi) ~ 3(h 2 ) 2 /2N, (h 5 ) ~ 2(h 2 )(h 3 )/5N, and 
(hg) ~ 5(n 3 ) 2 /289iV, respectively. Except for the nu- 
merical coefficient in the last relation, these results agree 
with the mean-field analysis in Ref. p0|| , although these 
approximate relations are not always satisfied by exact 
numerical results. 



V. CONCLUSIONS 

In this paper we have proposed an operator-algebraic 
approach to analyzing the yrast spectrum of weakly in- 
teracting bosons. The basic idea is to exploit the macro- 
scopic occupation of a particular mode to rewrite the 
second-quantized interaction Hamiltonian under the re- 
striction of Lhe lolal pui'liclo-immber conservation and 
the migulai-iiiuiiiL'iiliuii cunsei valiuii. 



We have analytically shown that the lowest-lying 
excitation spectrum from the yrast line is given by 
27(7713(71.3 — l)/34. Equation ( ^p|) also shows that 
this quasi-degenerate energy spectrum is separated from 
other levels by energy gaps that are linear in N. These 
results are in agreement with what we have found by 
numerical diagonalization of the same Hamiltonian [[16| . 
Equation (^) can also explain the excitations whose en- 
ergies are linear in N, again in agreement with our exact 
diagonalization calculations pl| . 



The method presented herein only exploits the pres- 
ence of a macroscopically occupied state, and it should 
therefore not be restricted to the case of L/N <C 1 dis- 
cussed above. Such an extension is in progress and will 
be reported elsewhere. 
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J L 


7 


7 


7 


8 


8 


8 


8 


excitation energy 


2502.8 


5000.0 


6253.3 


2507.7 


5000.0 


6254.0 


10628 



TABLE I. Excitation energies measured from the yrast line with N = 10000 for L — 7 and 8. The numerical results agree 
with our analytical predictions of N/4 = 2500, N/2 = 5000, 5iV/8 = 6250, and 17JV/16 = 10625 with O(iV ) deviations that 
arise from the coupling to the octupole mode through Eq. (ff5[). 
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